Percolation in a Class of Band Structured Random Matrices 
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We define a class of random matrix ensembles that pertain to random looped polymers. Such 
random looped polymers are a possible model for bio-polymers such as chromatin in the cell nucleus. 
It is shown that the distribution of the largest eigenvalue Amax depends on a percolation transition 
in the entries of the random matrices. Below the percolation threshold the distribution is multi- 
peaked and changes above the threshold to the Tracy- Widom distribution. We also show that the 
distribution of the eigenvalues is neither of the Wigner form nor gaussian. 
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Statistical properties of complex systems can be de- 
scribed by random matrix ensembles [l|. Most of the 
work has been concentrated on describing physical sys- 
tems where the resulting random matrix elements all are 
independent and identically distributed (iid) or ensem- 
bles of matrices where the probability distribution is in- 
variant with respect to orthogonal or unitary transforma- 
tions. These two possible generalizations of the gaussian 
unitary (orthogonal) ensemble represent different depen- 
dencies of the entries. While the first possibility has ar- 
bitrarily distributed entries, the second imposes strong 
statistical dependence over long distances resulting in dif- 
ferent properties of the eigenvalue spectrum 0, [1] . Com- 
mon to the random matrix ensembles based on the first 
possibility is the Wigner semicircle Q distribution for 
the eigenvalue spectrum and the Tracy- Widom distribu- 
tion [J] for the largest eigenvalue, even in the case of a 
band structured matrix [5|. 

For biological systems it turns out that random matri- 
ces are needed where all of the matrix elements are iid 




FIG. 1: Shown is the construction and a sample of the random 
matrices 



except for the main diagonal and first off-diagonal [6|, 
to describe, for example the conformations and three- 
dimensional organization of chromatin in the cell. In 
some case it is necessary to eliminate some matrix el- 
ements resulting in a band structured random matrix. 
Band structure random matrices have also been con- 
sidered for discretized models of solid state physics 0. 
Thus, while most entries are independent, the diagonal 
introduces correlations. It is thus natural to suspect that 
the resulting eigenvalue spectrum may not belong to ei- 
ther of the two universality classes. 

In this work we define an ensemble of band structured 
random matrices whose matrix elements can take on two 
possible values with the main diagonal matrix elements 
as dependent random variables. We are interested in the 
resulting eigenvalue spectrum and the distribution of the 
largest eigenvalue in terms of a percolation transition. 

Let h denote the bandwidth of a, N x N matrix with 
a displacement d from the main diagonal {N > 6 > 1) 
with matrix elements (see Figure [1]) 



(1) 



where 
Hjv(i,j) = for 



\i-j-d\>b 

\i- j\ <d,i^j,j^i-l,i + l 



and, within the band, we have Bernoulli variables 

h{i,j) = —1 with probability p £ [0, 1] 
j) — with probability 1 — p 



with the additional constraints 
h(i — 1, z) = h(i, « + 1) = — 1, 



2<i< N -1 



(2) 
(3) 

(4) 



N~"- is a norming factor for the eigenvalues. Thus the 
matrices that we consider are real symmetric band struc- 
tured random matrices. The band consists of random 
variables that take on the values —1,0 and the band may 
have been shifted. The main diagonal is the sum of the 
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FIG. 2: Finite-size scaling of the percolation probability for 
the ratio of matrix size A'^ and band width b of /N = 1. The 
exact values of the critical exponents for the two-dimensional 
site percolation were used. The non-universal value for pc 
was determined to be 0.61. Other ratios b'^ /N = const also 
produced good scaling. The data presented in the figure were 
obtained by averaging over at least 10000 random matrix con- 
figurations for each p and matrix size and band width combi- 
nations. 



off-diagonal random values with the exception of the first 
off-diagonal which is never shifted and has always values 
— 1. Thus the matrix is diagonally dominant. 

This structure mimics a polymer chain, where the first 
off-diagonal ensures the integrity of the chain, the diago- 
nal gives the coupling strength and the other off-diagonal 
elements mimic random loops. Note that the regular 
polymer case of 6 = 1, d = with gaussian random vari- 
ables has been solved by Dyson 3. 

Before we analyze the eigenvalue spectrum of the above 
defined random matrix ensemble, we look at the struc- 
tural properties of the random matrix. Since the ele- 
ments of the band are variables which can be with — 1 or 
it is tempting to suspect that the geometric distribu- 
tion of these has an influence on the eigenvalues. If we 
view the random matrix with the entries —1 or as an 
adjacency matrix representing an undirected graph (con- 
sidering only the upper triangular matrix, neglecting the 
diagonal and the first off-diagonal), then the graph may 
be a set of unconnected subgraphs or may have a per- 
colation structure. If the graph includes a percolating 
graph then this introduces a divergent length scale as 
N 00 possibly influencing the eigenvalues through the 
correlations introduced into the main diagonal. 

We define percolation in terms of a site percolation 
problem with free boundaries except for the left edge 
entries which are always fixed to —1 (see above). Here, 
the site percolation probability Poo is defined by the mass 
of the largest cluster divided by the size of the band. Due 
to the symmetry, we only consider the upper triangular 
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FIG. 3: Shown is the distribution of the eigenvalues. Note 
that we have shifted the distributions by their mean values. 
Indicated is also the Wigner distribution and a fit by a Gaus- 
sian. 



matrix (band). Two matrix elements having both a value 
of —1 belong to the same cluster if their indices differ at 
most by 1, i.e., are nearest neighbors in the sense of a 



simple square lattice [11 



To analyze percolation in the upper band, we use a 
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finite-size scaling ansatz 



(5) 



The finite size sca ling of the percolation probability (or- 
der parameter) [ij for this staircased slab with aspect 
ratio b'^ /N = 1 is shown in Figure[21 The critical value for 
the percolation threshold Pc for the simple square lattice 
is Pc = 0.592746 [lH and the exact values for the critical 
exponents are a = — 2/3, /3 = 5/36, = 4/3 As- 
suming universality we use these values for the critical 
exponents and find a value Pc — 0.61 for the percolation 
threshold that gives the best fit to a scaling function V. 
Of course, we don't expect Pc to be a universal quantity. 
Beside possible corrections to the finite-size scaling, the 
percolation problem posed by the staircased slab geom- 
etry is seen to be in the same universality class as the 
two-dimensional site percolation. 

Thus, at Pc we have a correlation length which is di- 
verging as the linear dimension of the matrix goes to 
infinity. It should be noted that if b ^ N, then many 
percolating clusters in the b direction can occur In 
the geometry considered here only one percolating cluster 
exists [l3| . In the following we explore the consequences 
of this on the eigenvalue spectrum and in particular on 
the distribution of the largest eigenvalue. 

Before we turn our attention on the distribution of 
the largest eigenvalue, we examine briefly the distribu- 
tion of the eigenvalues themselves. As suggested by the 
above percolation analysis we expect the geometric per- 
colation to influence the eigenvalue distribution. The 
distribution is shown in Figure [3] for three values of p. 
The data shown in the figures were obtained by averag- 
ing over 10000 Monte Carlo samples. Beside the matrix 
sizes shown, data were also produced for matrices rang- 
ing from N = 196 up to TV = 4225 for P/N = 1,5,10 
and d = 0, 10. 

First we note that for different combinations of N and b 
at fixed ratio b'^/N = 1 the data scales, in agreement with 
other band random matrix ensembles [5| . As a function of 
p the distribution undergoes a shape transition. The sign 
of the asymmetry changes and the distribution is nearly 
symmetric close to the percolation threshold. There we 
also compare the distribution to the Wigner distribution 
and to a Gaussian. Neither fits the distribution. Note 
that for the gaussian unitary ensemble a rapid change has 
been observed [1, , as one increases the band width, 
from a gaussian to a semicircle distribution. 

Let Amax(^^Ar) denote the largest eigenvalue of the 
random matrix Hn and 



'N,b,d. 



Pi 



N,b,d,p 



(A < Amax) 



(6) 



the corresponding distribution function. For the non- 
banded ensembles GUE and GOE, the distribution func- 
tion is the Tracy- Widom distribution 0]. This distri- 
bution is expected to hold for an even broader class of 



random matrices and the precise characterization of the 
class is now being investigated 15, [l3|. Known is that 
for symmetric matrices with iid entries of variance 1/A'^, 
such that all moments are finite, the Tracy- Widom dis- 
tribution holds asymptotically. If the entries decay with 
a power law then it is expected that the matrices belong 
to a different universality class 17[. 

To discover the phase transition in the behavior of the 
largest eigenvalue we calculate the cumulant of the dis- 
tribution of the largest eigenvalue and analyze the flow of 
the cumulant as a function of the matrix size at fixed ra- 
tio b'^ /N = 1 and the probability p with which the matrix 
elements are set. We define the fourth order cumulant as 



uiN,p)r^l-{X^)/mxY). 



(7) 



For a second-order phase transition the cumulant flows 
to three fix points [10|. Two fix points characterize the 
values of p below and above Pc and one characterizes the 
value of Pc- 

In Figure |4] is shown the fourth-order cumulant for the 
distribution of the largest eigenvalue as well as examples 
of the distribution above and below pc- While the distri- 
bution of the largest eigenvalue below Pc shows multiple 
peaks, above Pc these are reduced to a single peak. This 
change is also seen in the flow of the cumulant. For val- 
ues p well above the percolation transition point Pc the 
cumulant is almost constant and nearly zero. Below the 
percolation transition point the cumulant slightly fluctu- 
ates in a narrow band. Close to the transition point the 
flow changes strongly indicating the transition. 

In Figure |4] is also shown a fit of the Tracy- Widom 
distribution to the data above Pc at 0.9. The fit is very 
good indicating that in this ensemble above Pc the Tracy- 
Widom result holds. So far, for symmetric N x N ma- 
trices with iid entries of variance 1 /N, the Tracy- Widom 
distribution has been shown to hold [15| , while if the en- 
tries fall off with a power law, the matrices fall in a in 
different universality class [l9l |. Thus it is quite surpris- 
ing that for the class of random matrices defined here we 
do get a Tracy- Widom distribution. 

The geometric structure of a random matrix deter- 
mines the eigenvalues of a matrix to a large extend. In 
this paper we have introduced a class of band struc- 
ture random matrices that pertain to the description of 
polymers and may be capable to describe the observed 
three-dimensional organization of chromatin. Due to the 
Bernoulli type of its off-diagonal matrix elements we are 
able to define a percolation problem and link the change 
in the eigenvalues to the occurrence of a percolation tran- 
sition. The change is most noticeable in the density dis- 
tribution of the largest eigenvalue. The multiple peak 
distribution gives way to a single peak at the percola- 
tion transition. But also the eigenvalue density shows a 
change in skew at the transition point. Unclear is yet the 
type of distribution for the eigenvalues as well as for the 
largest eigenvalue below the transition. 



4 





0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 0.0045 0.005 0.0055 
N 




160 
140 
120 
100 
80 
60 
40 
20 





FIG. 4: Shown is the cumulant for the probabihty distribu- 
tion of the largest eigenvalue. The figure shows the cumulant 
UN,p = 1 - ((A - {A))*)/((A - {X)ff keeping the ratio b'^/N 
fixed to 1. Above the percolation transition at Pc (roughly 
at 0.61) the flow changes and for higher values of p is almost 
constant and zero. For values below pc the cumulant is fluc- 
tuating in a narrow band for the matrix sizes considered in 
this work. Also shown is the distribution at p = 0.1 and at 
p = 0.9. While below pc the distribution is multi-peaked we 
find a transition to a Tracy- Widom distribution above Pc- 
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